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DYNAMIC  STRESSES  IN  AN  ELASTIC 


CYLINDRICAL  LINING  OF  ARBITRARY  THICKNESS 


IN  AN  ELASTIC  MEDIUM 


ABSTRACT 


Dynamic  stresses  in  a  thick-wall  elastic  cylinder  in  an  infinite  elastic  medium 
during  passage  of  plane,  compressional  waves  are  investigated.  Dynamic 
stresses  around  the  cylinder  in  the  elastic  medium  are  also  determined.  Nu¬ 
merical  results  for  two  different  cylinders  with  ratios  of  outer  radius  to  inner 
radius  ranging  from  1.  05  to  1.  20  are  presented  in  a  dimensionless  form.  It 
is  shown  that  increasing  thickness  does  not,  in  general,  reduce  stresses  in 
the  cylinder;  in  addition,  dynamic  stresses  at  certain  wave  numbers  are 
higher  than  the  corresponding  static  value. 
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SECTION  I 


INTRODUCTION 


The  problem  of  dynamic  stresses  in  an  infinite  medium  containing 

cavities,  rigid  inclusions,  and  elastic  inclusions  has  been  studied  exten- 
[  1  2  3  4] 

sively  ’  ’  ’  the  problem  of  scattering  or  diffraction  of  sound  or  stress 

waves  by  thin  elastic  shells  in  fluid  and  elastic  media  has  also  been  the 

....  .  ..  ..  [5,  6,  7] 

subject  of  many  investigations. 


The  problem  presented  herein  is  the  response  of  a  thick-wall  cylinder 
in  an  infinite  elastic  medium  subjected  to  a  progressing  compressional  wave. 
It  is  assumed  that  the  thick-wall  cylinder  is  of  infinite  extent  and  is  embedded 
in  and  bound  to  an  infinite  elastic  medium.  A  plane  compressional  wave  of 
harmonic  time  variation  propagates  in  the  positive  X-direction  (Fig.  1 )  and 
impinges  on  the  embedded  cylinder. 


It  is  known  that  if  the  material  constants  (X,  p,  and  p)  of  the  cylinder 
and  the  elastic  medium  are  different,  scattering  will  occur  at  the  interface 
of  the  cylinder  and  the  surrounding  medium.  If,  however,  the  material  con¬ 
stants  are  the  same,  the  problem  should  reduce  to  the  simple  cavity  case 
presented  in  reference  1. 

In  general,  at  the  boundary  of  the  cylinder  and  the  surrounding  medium, 
the  incident  compressional  waves  are  reflected  and  refracted  as  compressional 
and  shear  waves.  Hence,  there  will  be  two  refracted  waves  propagating  into 
the  cylinder.  These  waves  will  be  reflected  at  the  traction-free  boundary  of 
the  inner  surface  of  the  cylinder.  Therefore,  seven  waves  exist;  one  incident 


1 


wave,  two  reflected  waves  in  the  elastic  medium,  and  four  in  the  cylinder  due 
to  refraction  and  reflection. 


The  solution  of  this  problem  involves  finding  the  coefficients  associated 
with  the  six  unknown  waves.  This  is  accomplished  by  four  equations  of 
continuity  at  the  interface  and  two  boundary  conditions  at  the  traction-free 
inner  surface  of  the  cylinder. 

!  !  I  i  i  I  1 


sFig.  1  Cylindrical  Lining 
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SECTION  II 


GENERAL  THEORY 


If  it  is  assumed  that  the  cylinder  is  infinite  in  extent,  the  problem 
becomes  one  of  generalized  plane  strain.  The  displacement  equation  of  motion 
is 


(X  +  p)vV 


u  +  p  v  u  =  p  u 


(1) 


where 

u  is  the  displacement  vector 
V  is  the  gradient  operator 
X  and  p  are  the  Lame  Constants 
p  is  the  density. 


The  displacement  vector  u  can  be  represented  by  a  scalar  potential 
and  a  vector  potential;  in  the  case  of  plane  strain  this  is 

u  =  V  <P  +  V  x(e  ip)  (2) 

z 


Each  potential  function  then  satisfies  a  scalar  wave -equation 

c  2  V  2  <P  =  '<P 
a 

2  2 

C/JV  $  =  ip 


(3) 

(4) 
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In  Equations  (2),  (3),  and  (4),  e  is  a  unit  vector  along  the  axis  of  the 

z 

cylinder  and 

2  X  +  2  u  2  _  M 

C  -  C  rt 

a  p  P  P 

As  shown  in  Fig.  1,  in  plane  polar  coordinates  (r,  6)  the  scalar  form  of 
Equation  (2)  is 


_  3  <p  1  3  £ 

Ur  3  r  +  r  3  6 


_  1_  3  0  _  3 

U0  r  3  9  3  r 


(5) 


and  the  stresses  are  related  to  the  potentials  by 


T  =  Av' 

rr 


+  i  ifjt. 

W  r  3r83  r2  3<7 


T00 


7  j. 

\r2  as2  r  3r  r  3r39  -2  86 


T  = 

r0 


o  u  /  i  a2»  i»J,i  »i.  i  sfi  l  „2  ... 

2Ml  r  Srde  2  sF  +  7  3t  +  1  ~2  '  J  V  *)  (6) 

'  r  r  d  0  / 
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SECTION  m 


INCIDENT  REFLECTED  AND  REFRACTED  WAVES 


For  convenience,  the  surrounding  elastic  medium  and  the  thick-wall 
cylinder  will  be  denoted  as  regions  Nos.  1  and  2,  respectively. 

The  incident  wave,  propagating  in  the  positive  X-direction  is  represented 

by 


<P 

$ 


(i) 

(1) 

(i) 

(1) 


i(Q'ix-  wt) 


0 


(7) 


where 

v 

is  a  measure  of  the  amplitude 
ui  the  circular  frequency 

a  =  lo  /  is  the  wave  number  of  the  compressional  wave. 
1  c 


In  polar  coordinates 


e  in  J  (ce ,  r)  cosn0  e 
n  n  1 


-iwt 


n=0 


(8) 
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where 


J  denotes  the  Bessel  function  of  the  first  kind  of  order  n 
n 


e  is  a  constant  so  that 
n 


1  n  =  0 


2  n  -  1 


Let  the  origin  of  the  polar  coordinates  coincide  with  the  central  axis  of  the 
cylinder.  The  waves  in  regions  Nos.  1  and  2  can  then  be  expressed  as  follows. 


REGION  NO.  1 


=  y  A  r)cosn0e  (9) 

(1)  Lj  n  n  i 

n  =  0 


*'(f)  =  I  BnH(n1)(g1r)sinnee-"4’t 

n  =  0 


(10) 


6 


REGION  NO.  2 


<P 


(1) 

(2) 


-l 


M  (a  r)  cosn0  e 
n  n  2 


n  =  0 


(1) 

(2) 


-l 


(2)  „  -iwt 

N  IT  '  (Br)  sinn0  e 
n  n  yp2  ' 


n=  0 


0 


(R) 

(2) 


■l 


R  (aor)cosn0e  iWt 
n  n  v  2  ’ 


n  =  0 


P 

\2) 


■I 


m  -ia>t 

S  H'  '  (/3  r)  sinn0  e 
n  n  2 


n  =  0 


(11) 


(12) 


(13) 


(14) 


In  Equations  (9)  through  (14) 


/T3\  /r>\ 

^(1)  ^  (1)  rePresent;s  compressional  and  shear  waves  in  region  No.  1 

^/o\  $  /o\  rePresents  the  inward  propagating  compressional  and  shear 

(2)  (2) 

waves  in  region  No.  2 

^(2)^  ^  (2)^  are  outSoinS  waves  in  region  No.  2 


A  ,  B  ,  M  ,  N  ,  R  ,  and  S  are  expansion  coefficients  to  be  determined 
n  n  n  n  n  n 


a  and  oi  are  the  compressional  wave  numbers  in  regions  Nos.  1  and  2 

J.  ^  % 
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P  and  B  are  the  shear  wave  numbers  in  regions  Nos.  1  and  2  (/3  =  —  etc.) 

1  c<3 

P1 

^  and  denote  the  Hankel  functions  of  the  first  and  second  kind  of 
n  n 

order  n  .  The  Hankel  function  of  the  first  kind  is  used  for  diverging  waves; 
e.g. ,  in  Equations  (9),  (10),  (13),  and  (14).  The  Hankel  function  of  the 
second  kind  represents  converging  waves;  e.g. ,  in  Equations  (11)  and  (12). 


In  Equations  (9)  through  (14)  there  are  six  undetermined  coefficients.  The 
boundary  conditions  which  allow  these  coefficients  to  be  determined  are  listed 
below. 

At  r  =  b  the  condition  of  continuity  requires  that 


T  =  T 

rr  rr 

(2)  (1) 


Tr  0  Trd 

(2)  (1) 


u  =  u 
r  r 

(2)  (1) 


U0  U0 

(2)  (1) 


(15) 


At  r  =  a  the  traction-free  boundary  implies 


T 

rr 


(2) 


=  0 


(16) 
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SECTION  IV 


SOLUTIONS 


isre  t  t  u  u  are  the  stresses  and  dis- 

rr(l),(2)  r  (1) ,  (2)  r(l),(2)  (1),(2) 

placements  due  to  the  total  displacement  potential  in  regions  Nos.  1  and  2;  e.  g. , 


T 

rr 


(2) 


(17) 


where 


cf)  =  <pW  + 

(1)  (1)  (1) 

*«-«&? 


<p  =  <f>W  +  cp^ 
(2)  r(2)  y(2) 

*  (2)  *  (2)  +  *  (2) 


(18) 


Substituting  Equations  (9)  through  (14)  and  (18)  into  Equations  (5)  and 

(6)  yields  the  corresponding  displacement  and  stress  components  in  regions 

— i  co  t 

Nos.  1  and  2.  With  the  time  factor  e  omitted,  the  expressions  for  the 
displacement  are 


9 


u  =  r 

r<D 


0° 

n=0 


cosnfl 


u  =-r 

(1) 


o° 

1  I  [Vn‘nnVV)+ VHn>(V)  + Vl^n’  ^  _ 


n  =  0 


(19) 
sin  n0 

(20) 


°o  r 

ur  =r'1  I  LMnVH<n'<a2r>  +  NnnHn)(V)  +  RnVHn)  (“2r> 

(2)  n=0 


r(l) 


+  S  nH  (j8  r)  cos  nfl 
n  n  2  1 


] 


(21) 


u. 


=  -r 


-1 


(2) 


oo 

n 

n=  0 


M  nH^  (a  r)  +  N  /3.rH^  (lr)  +  R  nH^V0r) 
n  n  2  '  n2  n  2  '  n  n  2  ' 


+  S  firH(1)'  (13  r) 
n  2  n  2  J 


sin  n0 


(22) 


The  expressions  for  the  stress  are 


T 

rr 


(1) 


=  2Mlr 


OO 

I 

n  =  0 


(4>  e  in  ,DW  +  A  D'*v  +  B  . 
On  lnr  -  n  1  nr  nl 


D(R)  +  B  $  (R* )  cos  n 6 
nr  nl  nr 


(23) 


an  =  2**ir 

9V) 


-2 


(tf>  e  in.F 

'On  1  nr 


<‘)  + 


F<R>  -  B 


nl  nr 


$(R) 


nl  nr 


)  cosnfl 


(24) 


n=0 
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rO 


(1) 


oo 

=  2M  r"2  y  (4>  e  inR^+A  ,E(R)  +  B  &(R))sinn0  (25) 

1  Lj  On  1  nr  nl  nr  nl  nr' 


n=  0 


00 


rr 


(2) 


2V 


-2 


(M  0D(1)+N  +  R  n<R>  +  S  o^(R))cosn0  (26) 

n2  nr  n  2  nr  n  2  nr  n2  nr  '  ' 


n  =  0 


00 


(2) 


w 

=  2M0r"2^)  (M  „F(1)-N  $(l)  +  R  F(R)  -  S  $(R))cosn0  (27) 

2  L,  n2  nr  n  2  nr  n  2  nr  n  2  nr  '  v  ’ 


n=  0 


r  0 


(2) 


2V 


oo 

y  (M  E(l)  +  N  6(l)  +R  0E(R)  +  S  og(R))sinn0  (28) 
Lj  n  2  nr  n  2  nr  n  2  nr  n  2  nr 


n=  0 


with 


,D^  =  (n2  +  n-^-  ^  r2)  J  (a  r)  -  ot  r  J  .(err) 
1  nr  2  1  '  n  1  '  ln-11' 


(i)  2 

,EW  =  (n  +n)J  (a ,r)  -narj  .(err) 
1  nr  n'  1  1  n-1  1 


lFnr  "(n2  +  «  -  “i  *  +  \  f<\  rVn(  V>  +  VJn-l<V> 


(29) 
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.D(R)  =  (n2  +  n  -  ^  2  r2)  H^(a.r)  -  a.rH^  (a.r) 
3  nr  '  2  j  '  n  '  3  '  j  n-r  3  ' 


!<*>  = 


r(!) 


.E  =  (n  +n)H  (a.r)  -  nff.r  H  ,  (a.r) 
3  nr  n  ]  3  n-1  3 


F(R) 
j  nr 


2  2  2  1  2  2  fl\  n\ 

=  -  (n  +  n  -  a  .  r  +r  r  )Hu(ff.r)  +  a.rHu1(a.r)  (30) 
3  2  3  n  3  3  n-lv  3 


j  =  1,  2 


(i)  2  12  2  (2) 

DW  =  (n  +n-— /3„r  )  H'  ' 


r(2) 


2  “nr  <n'+n'2  V  )Hn  (V>  “  VHn-l(V> 


E(l)  =  (n2  +  n)  H(2)  (a  r)  -  no;  r  H(2) ,  (a  r) 
2  nr  v  ’  n  v  2  '  2  n-lv  2  ' 


F^  = 
2  nr 


2  2  2  1  2  2  (2^  (2) 

=  -(n  +n-«2r  +  —  j3^r  )  IT  '  (a2r)  +  oi^r  _  1(a2r)  (31) 


(j)  (R)  _  _(n2  +  nj  H(1)  (p  +  np  r  H(!)  (p  rj 

3  nr  n  ^3  3  n-lv  3  ' 


&  ^  =  -(n2  +  n  -  ^  /3  2  r2)  (fi  .r)  +  /3.r  (0.r)  (32) 

nr  2  3  n  3  j  n-1  3  '  ' 


3  nr 

j 


=  1,  2 
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"  i^i  n  tn.n  •  fi^mpip  - 


2®ir  '  -<n2  +  n>  Hf  <V>  +  nV  Hn- l^jr> 


!®nr  '  -<n2+n‘|'3lr2)Hi2>^2r)  +  ^2rHi2-l^2r>  (33> 


and 


dH 


H'  (x)  =  ~~ 
n  dx 


dJ 

J'  (x)  =  . 

n  dx 


n 


(34) 


Equations  (19)  through  (28)  are  the  general  expression  for  the  stresses 
and  displacements  in  regions  Nos.  land  2.  To  evaluate  the  coefficient  A  ,  etc., 
use  the  conditions  of  continuity  at  r  =  b  and  the  boundary  conditions  at  r  =  a  . 

At  r  =  b 


T  =  T 

rr  rr 

(2)  (1) 


M  D(l>  +  N  $(l>  +  R  9D<R>+S  $(R)  -  v(A  D(R)  +  B  <t>  e  in  D 

n2  nb  n2  nb  n2  nb  n2  nb  nl  nb  nl  nb  On  1  nb 

(35) 


ry(2)  r0(l) 


M  E(l^+  N  +  R  ,E<R>  +  S  g(R)  -  u(A  E(R)  +  B  6<R))  =  v  cp  e  i"  E 

n2nb  n2nb  n2nb  n2  nb  '  nl  nb  nl  nb  On  1 


{n  r(i) 
nb 


(36) 
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u  =  u 
r  r 

(2)  (1) 


M  ttbH(2}'(ab)  +  N  n  H(2)(d  b)  +  R  a  bH(1)'(ab)  +  S  nH(1)  (fib) 
n  2  n  2  n  n  2  n2n  2  nn  2 

-  [  AnaibH^1)’  (o^b)  +  BnnHji1)  (^b)  ]  =  ^i"  a^bJ^b)  (37) 

U0  =U0 

(2)  (1) 

MnnHn  }  (a2b)  +  V2bHi2)'  (^2b>  +  VHn  *  (“2b)  +  WHn  ''  (/32b) 


-  |  A  nH(1)(a  b)  +  B  ^bH(1)'(Sb)  1  = 
Ln  n  1  n  1  n  '  1  '  J 


</>  e  innJ  (0!  b) 
On  n'  1  ' 


(38) 


At  r  =  a 


rr 


=  0 


(2) 


M(2>  +  N  +  R  9D(R)  +S  ®(R) 

n  2  na  n  2  na  n  2  na  n  2  na 


=  0 


(39) 


r0 


=  0 


(2) 


M  E*1*  +  N  S(l)  +  R  +  S  S(R) 

n  2  na  n  2  na  n  2  na  n  2  na 


=  0 


(40) 


where  v-  u  /u  is  the  ratio  of  the  shear  moduli  of  regions  Nos.  land 2. 

-L  Z 


Thus,  the  coefficients  A  ,  B  ,  M  ,  N  ,  R  ,  and  S  are  determined  by  the 

n  n  n  n  n  n 

six  simultaneous  Equations  (35)  through  (40). 
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Using  matrix  notation,  Equations  (35)  through  (40)  can  be  conveniently 
expressed  as  follows 

MW  -w 

where 


,d<r> 

1  nb 

1  nb 

D(i) 

2  nb 

2)(i) 

2  nb 

d(R) 

2  nb 

$(R> 

2  nb 

iR) 
v.  Ev  , ' 

1  nb 

,g(R) 

1  nb 

E(i) 

2Enb 

c(i) 

2  nb 

f(r) 

2  nb 

c(R) 

2  nb 

QbH(1)(ffb) 

1  n  1 

nH^flJjb) 

“2bHn2)'“2b) 

nHn^2b> 

“2bHn1)<“2b) 

"Hl1>^2b> 

nH(1\ab) 
n  1 

^1  bHn)(^ib) 

“Hn2,<“2b) 

'32bHnWb> 

nH^V  b) 

n  z 

^2bH<„l>^2 

0 

0 

D(l) 

2  na 

2  na 

D(r> 

2  na 

®<R> 

2  na 

0 

0 

,E<‘> 

2  na 

«« 

2  na 

e(r) 

2  na 

26<R) 

2  na 

(42) 
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and 


Hence 


(43) 


(44) 


(45) 


where  [A„]  is  the  adjoint  of  [a_]  and  |a|  is  the  determinant  of 


16 


t 


REDUCTION  TO  THE  SIMPLE  CAVITY  CASE 


Consider  the  case  when  the  material  constants  of  the  two  regions  are  the 


same: 


Therefore 


It  follows  that 


WA 


"i =  ^  “ 


p  =  p  =  p 
H1  2 


a  =  a  =  cl 
12 


^1  ^  2  ^ 


^  =  1 


d(R)  = 

1  nr 


d(R) 

2  nr 


$<R>  = 
1  nr 


,®(R) 

2  nr 


lE<R>  - 

1  nr 


yR> 

2  nr 


ig 

1  nr 


(R)  _ 


g(R) 

2  nr 


(46) 
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Furthermore,  using  the  relationship  between  and 

n  n 


H(1)(x)  +  H(2)(x)  =  2J  (x) 
n  n  n 


(47) 


it  can  be  shown  that 


9d«  +  9d<r> 

2  nr  2  nr 


=  2  D(l) 
1  nr 


(48) 


E(i)  +  e(R)  =  2  E(i) 
2  nr  2  nr  1  nr 


(49) 


Substituting  Equations  (46),  (48),  and  (49)  into  Equation  (45),  the 
undetermined  coefficients  for  the  two  potentials  in  region  No.  1  are 


A 

n 


-  4>  e 
0  n 


.n 

l 


D(i) 

3)<r> 

1  na 

1  na 

E(i) 

g(R) 

1  na 

1  na 

d<R) 

®(R> 

1  na 

1  na 

e(r) 

g(R) 

1  na 

1  na 

(50) 
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and 


,d(R) 

1  na 

1  na 

<R> 

E^l) 

1  na 

1  na 

d(R) 

®<R) 

1  na 

1  na 

,e<r> 

,s<R> 

1  na 

1  na 

(51) 


These  expressions  for  A  and  B  are  the  same  as  those  in  reference  1.  The 

n  n 

remaining  coefficients  M  ,  N  ,  R  ,  and  S  are 

n  n  n  n 


M 

n 

N 

n 


1 

2 

0 


.n 

i 


S 

n 

R 

n 


B 

n 


1  a 

- ■  <p  e  i 

2  r0  n 


n 


+  A 


n 


(52) 


Thus  the  shear  potential  in  region  No.  2  is  obviously  the  same  as  that  in 
region  No.  1.  The  total  scalar  potential  in  region  No.  2  is 
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*<2>  “  ^(2)  +  *<2)  -  I  2  Vn1”  H<n2>  (“r)oos  n9 


n=0 


(  b  0n  e  in  +  A  )  (ar)  cos  nfl 
2  0  n  n  n 


OO 

'  1  kvnj„<“r>  + 

n  =  0 


cos  n  9 


which  is  identical  to  the  scalar  potential  in  region  No.  1.  Therefore,  the 
problem  is  reduced  to  a  simple  cavity  case. 


(53) 
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SECTION  V 


NUMERICAL  RESULTS  AND  DISCUSSIONS 


Numerical  results  are  obtained  for  r  and  r  at  r  =  b  and  t 

i*i*  w  u  u  u 

(1)  (1)  (2) 

at  r  =  a.  This  is  accomplished  by  summing  the  series  in  Equations  (23),  (24), 
and  (27),  respectively. 

For  reasons  of  expediency,  the  expressions  for  the  stresses  are  non- 
dimensionalized.  The  important  nondimensional  parameters  can  then  be 
properly  identified,  and  the  effects  of  these  parameters  on  the  stresses  in  the 
cylinder  as  well  as  in  the  surrounding  medium  can  be  evaluated. 

To  nondimensionalize  these  expressions,  recall  the  expression  for  the 
incident  wave  given  in  Equation  (7)  as 


V 


i(aix  -  wt) 


Therefore,  the  stress  in  the  incident  wave  in  the  direction  of  propagation  is 


T 

XX 


2  -  wt) 

+  2MX)  4>0  e 


(54) 


2 

It  follows  that  -(A.^  +  2^)0^  is  the  stress  amplitude  of  the  incident  wave. 
Denote  it  as 


T 

0 


+  2/^)  a 


Vo 


(55) 
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Thus,  Equations  (23),  (24),  and  (27)  can  be  nondimensionalized  by  dividing 
through  by  as  shown  below: 


rr 


111 _ 2. 


rr 


(1)  0 


<p  )32r2 
cri  n=  0 


5  (<V  in  D(l)  +  A  D1 

L  On  1  nr  nil 


(R) 

nr 


+  B  ®(R))cosne 
n  1  nr 


(56) 


T*  =  _!!ai  = 

00(1)"  to 


<£  /J2r2 
0P1  n  =  0 


y  (<P  e  i“  F(l)  +  A  F(R) 
Li  0  n  1  nr  n  1  nr 


-  B  )  cos  nd 

n  1  nr  ' 


(57) 


ee 


121 


ee  t 

(2)  o 


0  1  n  =  0 


y  (M  F(l)  -  N  ®(i) 
L  n2  nr  n  2  nr 


+  Rn2Fir)"S  o®(R))cosn0 
n  2  nr  n  2  nr  ' 


(58) 


where  v  is  the  ratio  of  shear  moduli  of  regions  Nos.  1  and  2. 


22 


By  letting  r  =  b  in  Equations  (56)  and  (57),  and  r  =  a  in  Equation  (58) 
we  have,  for  example, 


rr. 


(1) 

r  =  b 


vJb* 


2  (+„E  i"  •  D1'*  +  A 


On  1  nb 


nlD!b>+BnlIib>oosnl’  (59) 


n  =  0 


with  similar  expressions  for  Equations  (57)  and  (58) . 


By  examining  the  arguments  in  . . .  etc.  ,  that  are  associated  with 

various  terms  in  Equation  (59)  etc. ,  it  is  found  that  five  essential  nondimensional 
parameters  exist.  These  are 


v  = 


c 

a  a 
2  1 


a  c 

1  a. 


ratio  of  the  dilatational  phase 
velocities  of  regions  Nos.  1 
and  2 


b/ a 


ratio  of  dilatational  phase 
velocities  to  distortional  phase 
velocities  in  regions  Nos.  1 
and  2 


V 
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ratio  of  outer  radius  to  inner 
radius  of  the  cylinder 


By  substituting  these  expressions  into  Equation  (59)  . . .  etc. ,  the  final 

expressions  for  r  *  t*  at  r  =  b  and  r  *  at  r  =  a  are  obtained 
rr  99  99 

(1)  (1)  (2) 


T* 

rr 

(1) 
r  =  b 


,  , l  2  2,  v2  ^ 
n=0 


00 

N  (<t>  c  in,D^  +  A  D^  +  B  $(R))cosn0 
/_j  x  0  n  1  nb  nl  nb  n  1  nb 7 


(60) 


T* 

ee 

(i) 

r  =  b 


,  .  2  2,  .2  _ 

Vl  ’  <“ia)  n=0 


y  (4  £  i”,F*1*+A  F(^-B  2>^)cosne 

Lj  On  1  nb  n  1  nb  n  1  nb 


(61) 


"(2) 
r  =  a 


Vkl(“la>2  n=0 


Y  (M  FW  -N  +  R  F(R) 

Lj  n  2  na  n  2  na  n  2  na 


S  ®(R))cosn0 
n2  na' 


(62) 


The  expressions  for  can  be  expressed  in  terms  of  the 

1  nb  2  na 

incident  dimensionless  wave  number  (a^a)  and  the  parameters  defined  above. 
For  example 


lDib  -  O'2  +  "  -  \f>\  b2>  Jn(“1b)  ‘  “lb  Jn  -  l<“lb> 


2  12  2  2 

n  +  n  ~  2  kl  17  (aia)  J  Jn(T?Qfia)  “  ^“l^n  -  i(TJQ!ia) 


(63) 
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It  is  apparent  that  the  stresses  in  the  cylinder  as  well  as  in  the  surrounding 
elastic  medium  are  not  only  a  function  of  the  incident  wave  frequency,  but  also 
depend  on  the  four  physical  parameters  of  the  two  media  and  a  geometrical 
parameter  of  the  cylinder. 

By  restoring  e  -ia^,  Equations  (60)  through  (62)  have  the  form 


(R  +  il)  e 


-icot 


,r2  +  i2)  1/2  e -!(<*-«) 


(64) 


The  real  part  represents  the  stresses  at  t  =  0  and  the  imaginary  part 

gives  the  stresses  at  t  =  T/4,  where  T  =  2?r/w  is  the  period  of  the  incident 

2  2  l/2 

wave.  The  absolute  values  (R  +  I  )  correspond  to  the  maximum  values 

of  r*  ,  t*  and  t*  which  occur  in  the  interval  t  =  0  to  t  =  T/4  . 
rr  99  99 

(1)  (1)  (2) 

The  phase  angles  are  given  by  6  =  arctan  i/R . 


Numerical  results  are  obtained  for  two  cases  corresponding  to  a  soft 
and  a  stiff  cylindrical  lining. 

SOFT  CYLINDRICAL  LINING 

In  this  case  the  dimensionless  physical  parameters  are  assumed  to  be 


V  = 

2.90 

y  = 

1.50 

a  = 

1 

.25 

<7  = 

2  - 

.20 

T]  =  1.  05,  1. 1  and  1. 20 
The  results  are  shown  in  Figs.  2  through  6. 
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2 


Fig.  2  Distribution  of  Normalized  Tangential  Stress  at 
r  =  a  and  V  =  1 . 1  (Case  I ) 


Fig.  2  shows  the  angular  distribution  of  Jr*  j  for  77=1.1  at  r  =  a 

UU 

(2) 

for  two  wave  numbers:  a^a  =  .  20  and  a^a  =  1. 0.  At  the  low  wave  number, 
the  distribution  is  nearly  symmetrical  with  respect  to  the  y  axis;  at  o^a  =  1.  0, 

the  peak  stress  is  shifted  toward  the  incident  side  of  the  cylinder.  This  is  also 
found  in  reference  2. 


Fig.  3  shows  the  stresses  as  a  function  of  a^a  and  17  at  t  =  0  and  t  =  T/4. 
For  this  case,  an  increase  in  r\  does  not  change  the  stress  appreciably;  in 
fact,  as  rj  increases,  the  stress  also  increases. 


Fig.  4  shows  the  stresses  at  6  =  n/2,  n  and  r  =  a  for  .  10  ^  a  ^2.  0. 
The  maximum  stresses  at  7r/2  occur  at  a^a  «  .25;  this  agrees  with  the  results 
of  reference  1. 
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2.0 


0  1.0  20 


a, a 

Fig.  3  Normalized  Real  and  Imaginary  Tangential  Stresses  at 

7 r 

r  =  a,  0=-r—  for  Various  7]  (Case  I ) 


Figs.  5  and 6  show  the  effects  of  77  on  the  stresses  in  medium  (1)  at  r  =  b. 

It  is  seen  that  |r*  |  at  7r/2  decreases  as  tj  increases,  with  the  peak  value 

du 

(1)  I  I 

again  occurring  at  a  a  «  .25.  On  the  other  hand,  |  r*  J  in  all  cases  in- 

rr(l) 

creases  as  77  increases.  This  is  apparent  since  the  rigidity  of  the  cylinder 
is  increased  as  77  increases;  hence,  higher  radial  stresses  are  produced  at 
the  boundary. 
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0  1.0  2.0 

ata 


Fig.  4  Normalized  Tangential  Stress  at  r  =  a,  0  = 
for  Various  i)  (Case  I) 

STIFF  CYLINDRICAL  LINING 

The  dimensionless  parameters  are  assumed  to  be 

v  -  .  31 
y  =  .70 
<*1  =  .25 

a2  =  -3° 

tj  =  1. 05,  1. 1  and  1. 2 

The  results  for  this  case  (Case  II)  are  shown  in  Figs.  7  through  11. 
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4.0 


0  1.0  2.0 

a,a 


Fig.  6  Normalized  Radial  Stress  at  r  =  b,  0  = 


for  Various  rj  (Case  I ) 
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Fig.  7  shows  the  angular  distribution  of  |  r*  )  for  rj  =  1. 1  at  r  =  a 

(2) 

for  two  wave  numbers;  o^a  =  .20  and  a^a  =  1.0.  Note  that  the  shifting 
exhibited  in  the  case  of  the  soft  cylindrical  lining  (Case  I)  also  occurs  in  this 
case.  The  magnitude  of  the  stresses  are,  however,  much  higher. 

Fig.  9  shows  the  effects  of  rj  on  stresses  in  the  cylinder.  In  this  case, 
the  stress  decreases  as  rj  increases;  this  is  in  contrast  to  the  previous  soft 
cylinder  case. 


if 

Z 


Fig.  7  Distribution  of  Normalized  Tangential  Stress  at 
r=aandrj=  1.1  (Case  II) 
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10.0 


0  1.0  2.0 

a,a 


7T 

Fig.  8  Normalized  Real  and  Imaginary  Tangential  Stresses  at  r=  a,  9  =  — -  , 

t-d 

for  Various  tj  (Case  II ) 

Figs.  10  and  11  show  the  effects  of  r\  on  the  stress  in  region  No.  1.  The 
general  trend  is  the  same  as  in  (Case  I). 

The  contrast  in  the  results  for  the  stiff  lining  and  the  soft  lining  should 
be  emphasized.  Both  cases  exhibit  the  tendency  for  stress  to  concentrate  at 
large  wave  numbers  on  the  incident  side  of  the  cylinder;  however,  the  tangential 
stress  in  the  lining  is  significantly  higher  in  the  case  of  the  stiff  liner.  The 
wave  numbers  for  maximum  stress  at  9  =  7r/2  and  ir  are  approximately  the 
same  for  both  the  stiff  and  the  soft  lining,  but  the  effects  of  increased  lining 
thickness  on  tangential  stress  are  opposite  in  the  two  cases.  The  stress 
decreases  with  thickness  for  the  stiff  liner,  but  increases  with  thickness  in  the 
case  of  the  soft  liner. 
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0  1.0  2.0 

a,  a 


7 r 

Fig.  9  Maximum  Normalized  Tangential  Stress  at  r  =  a,  9  =  — —  , 
for  Various  n  (Case  II ) 


0  1.0  2.0 

a,  a 


Fig.  10  Normalized  Tangential  Stress  in  Medium  (1)  at 

7T 

r=  b ,6=  —  ,tt  for  Various  T)  (Case  n) 


7T 
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Other  points  of  contrast  are: 


(a)  The  variance  of  tangential  stress  in  the  cylinder  with  thickness 
is  less  in  the  case  of  the  soft  liner. 

(b)  The  tangential  stress  in  the  infinite  medium  at  the  boundary  of  the 
lining  is  less  in  the  case  of  the  stiff  lining.  In  addition,  increasing 
thickness  markedly  reduces  this  stress,  significantly  more  so  than 
for  the  soft  lining. 

(c)  The  radial  stresses  in  the  infinite  medium  at  the  boundary  of  the 
lining  are  much  less  in  the  case  of  the  soft  lining.  In  both  cases, 
this  stress  reduces  rapidly  with  decreasing  thickness. 


o  i .0  2.0 


Fig.  11  Normalized  Radial  Stress  in  Medium  (1)  at 

7T 

r=  b ,$= - ,7r  for  Various  rj  (Case  II) 

2 
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APPENDIX 


THE  EXACT  TWO-DIMENSIONAL  SOLUTION  FOR  AN 
ELASTIC  CYLINDRICAL  LINING  IN  AN  INFINITE  ELASTIC 
MEDIUM  UNDER  BIAXIAL  COMPRESSIVE  LOADING 

This  appendix  derives  expressions  for  the  stresses  in  an  elastic  cylin¬ 
drical  lining  in  an  infinite  elastic  medium  under  biaxial  compressive  forces. 
This  static  solution  represents  the  infinite  wave  length  (zero  frequency)  solution 
to  which  the  traveling  wave  solution  must  converge  in  the  limit.  This  is 
accomplished  by  using  a  superposition  approach  in  which  the  compatibility 
equations  for  the  stress  and  displacement  at  the  interface  between  the  lining 
and  the  infinite  medium  are  expressed  in  terms  of  unknown  radial  and  tangential 
tractions  assumed  to  act  on  the  lining  and  the  boundary  of  the  infinite  medium. 
The  solution  derived  is  a  plane  strain  solution. 

STATIC  STRESS  AND  DISPLACEMENT  OF  AN  UNLINED  CYLINDRICAL 
CAVITY  BOUNDARY  DUE  TO  THE  BIAXIAL  COMPRESSIVE  FIELD 

T  =  -T,T  =  e  T 

xx  0  yy  0 

Consider  an  infinite  elastic  medium  with  a  cylindrical  cavity  of  radius 

"b"  under  biaxial  compression  t  =  -  t  t  =  e  t  where  e  =  ~  cr  /l  -  or 

xx  0  yy  0  11 

(Fig.  12).  The  subscript  1  is  used  to  identify  the  infinite  medium  as  opposed 

to  the  lining  considered  below.  The  displacements  at  the  cavity  boundary  due 

[  7] 

to  the  biaxial  compressive  field  are 


t  b 

—  (l  +  o)  1  +  2(1  -  2<x  )cos  20 

E1  1  L  1 


(65) 
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1  .1 


l-o; 


Fig.  12  Unlined  Cavity 


u 


2r0b 


e. 


(i) 

r  =  b 


1  L 


(1  -  2 (Tj)  (1  +  a  )  sin  2 Q 


(66) 


Reference  8  gives  the  tangential  stress  at  the  boundary  of  the  cavity  due 
to  the  biaxial  compressive  field;  i.  e. , 


0 


ee ...  1-ct. 


-1  +  (2  -  4o^)  cos  2  9 


-(l)  "  l  L 

r  =  b 

The  radial  and  shear  stresses  vanish  at  the  boundary. 


(67) 
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STATIC  STRESS  AND  DISPLACEMENT  OF  AN  UNLINED  CAVITY  BOUNDARY 

DUE  TO  APPLIED  BOUNDARY  TRACTIONS  f  AND  r  n 

rr,,,  r  Q,,y 


Consider  an  infinite  elastic  medium  with  a  cylindrical  cavity  of  radius 

"b"  under  action  of  applied  radial  and  tangential  boundary  tractions  r  and 

rr(l) 

r  ,  respectively.  Let 
r  (1) 


Fig.  13  Boundary  Tractions  -  Unlined 
Cavity 
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X 


where  t  ,  t  ,  and  t  are  arbitrary  constants.  Fig.  13  illustrates 

rro  rri  1 

the  problem.  At  the  boundary  (r  =  b)  displacements  and  stresses  are 


r  =  b 


2(1  +  a  )b 


(1  +  a  )  + - 

1  V  E 


-  a  )  +  r 
*6  V  rfl 


cos  20 


(69) 


(1  +  c±)b 


sin  2  0 


(70) 


r  =  b 


and 


r  =  b 


T 

rro 


+ 


cos  2  0 


r  =  b 


-2  t  )  cos  2  0 
r0„  ' 


r  =  b 


sin  2  0 


(71) 
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STATIC  DISPLACEMENTS  AND  STRESSES  IN  AN  ELASTIC  CYLINDRICAL 
LINER  UNDER  ARBITRARY  BOUNDARY  TRACTIONS 


The  generalized  stress  function  and  stresses  for  plane  strain  in  cylindrical 
coordinates  (Fig.  14)  are 


^(2)(r»  0)  =  C1r2  fnr  +  C2r2  +  C^nr  +  +  (c^  +  Cgr4  +  — +  C^cos  2  6 


(72) 


C3  /  6C7  4C8  \ 

Trr  =  C1(l  +  2£nr)+2C2+—  -  1 2Cg  +  —  +— Icos20  (73) 
(2)  r  '  r  r  • 
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(2) 


°3  l  2  6C7  \ 

C  (3  +  2£nr)  +  2C2 - +  I  2Cg  +  12C6r  +  —  I  cos20 

r  \  r  > 


(74) 


r  6 


(2) 


2C  +  6C  r2 - 7 

v  5  6  r4 


6C„  2C, 


sin  2  6 


(75) 


where  the  subscript  2  identifies  the  elastic  lining  as  opposed  to  the  elastic 
infinite  medium. 
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Fig.  14  Elastic  Cylinder  with 
Boundary  Tractions 


The  stress-displacement  relationships  are 


du 


r  1  +  or 

J2I _ 2_ 


rr 


(2) 


3r 


(1  -  C T  )  T  -  <J  T 

2  rr(2)  2  W@>J 


u 


9u, 


=  _I2L  +I  _ (21  = 

00(2)  r  r  80 


1  +  g2 


(1-g2)T00  'Vrr 

(2)  rr(2)  J 


Integrating  these  equations  for  u  and  u  and  neglecting  rigid  body 

r(2)  (2) 

motions,  we  obtain 


(76) 
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1+<72^ 

ur  =  —  '|(1  "  2a2)  |  C1(2ri!nr  -  r)  +  2C„r 


(2)  2 


- 2a  C  r 

r  2  1 


C  2C 

-d-V4f-  +  4"2V3--5-  +  2C5r 


cos  2  0 


(77) 


1  +  a_ 


u. 


(2) 


j4(l-a2)Clr0+-r 


(12  -8a2)C6r 


+  4C. 


^  -  2<V  C8 

2 

r 


sin  2  0 


I 

\ 


(78) 


In  the  u 


necessary  that 


equations,  note  that  to  assure  periodicity  in  0,  it  is 


C  =  0  (79) 


APPLICATION  OF  BOUNDARY  AND  COMP  AT  ABILITY  EQUATIONS  TO 
DETERMINE  UNKNOWN  COEFFICIENTS 


The  system  of  equations  developed  above  contains  nine  unknowns;  i.  e. , 

C„,  C  .  C  ,  C_,  C  ,  t  ,  r  ,  and  r  .  Evaluation  of  these  con- 
2  3  5  6  1  8  rrQ  rr^  r01 

stants  entails  use  of  the  boundary  conditions  at  r  =  a  and  the  stress  and 

displacement  compatability  requirements  at  r  =  b . 
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Boundary  Conditions  at  r  =  a 


rr 

(2) 
r  =  a 


C  /  6C  4C  \ 

2C2  +"f  -  2C5  +  4  +  2  00829  =  0 

a  \  a  a  / 


6C 


r0 

(2) 

r  =  a 


2C  +  6C  a  - 
5  6 


a 


2C, 


a 


sin  2  0  =  0 


(80) 


Therefore 


2C2  4  -  =  0 
a 


(81) 


6C  4C . 

a  a 


(82) 


2  6C7  2C8 

2C5  +  6C6a  -  T  =  0 

a  a 


(83) 


Stress  Compatability  Relations  at  r  =  b 


rr 

(2) 

r  =  b 


c,  /  6C7  4Cs  \ 

2C  +  -  2C_  +  — 7  +  — —  cos  2  0  =  r  +  r  cos  2  0 

2  b2  V  5  b4  b2/  rro  rri 


6C, 


2C, 


r0 

(2) 

r  =  b 


2V6C6b  - 


sin  2  0  =  t  sin  2  0  (84) 

rSi 
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Therefore 


2C.  +  ~~  =  r 

2  b2  rro 


(85) 


6C_ 


-  2CL  + 


b2/  rrl 


(86) 


6C 


2C  +  6C  b  - 
5  6 


2C. 


8 

2  Tr0 


(87) 


Displacement  Compatability  Relations  at  r  =  b 


u  =  u  =  u '  +  u 

r  r  r  r 

(2)  (1)  (1)  (1) 

r  =  b  r  =  b  r  =  b  r  =  b 


u  =  u  =  u  ’ 
Q  Q 


+  u 

(2)  -(1)  V)  V) 


r  =  b  r  =  b  r  =  b  r  =  b 


(88) 


Using  Equations  (65),  (66),  (69),  (70),  (77),  and  (78),  there  results 


l  +  o\ 


<l-2a2)(2C2b)-  — 


rb 

X  <1  +  <v 


t  b 

rro 

-ir<1+(v 


(89) 
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a  +  <y 
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2  o 


2c 

+  2C5b- 


^8 
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2rob 


(l-2ff1)(l+ff1) 


2(1  -  oy 


(-  -  a  )  + 

v6  r 


r0 


(90) 


1  +  a2 
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(6  -  4a2)  C6b3  +  2C,_b  +  2 

b 


^  -  2VCg 

b 


2Tob  r 

—  [(1  -  2a1>(1  +  vj 


(1 +  <9 

E. 


2c r  \  +  r  (  — 
1;  r01l3 


-  2ax) 


(91) 


COMPUTATION  AND  NUMERICAL  RESULTS 

Equations  (81),  (82),  (83),  (85),  (86),  (87),  (89),  (90)  and  (91)  form  a 
system  of  nine  equations  for  the  nine  unknowns.  For  compatibility  with  the 
nondimensional  dynamic  analysis  presented  previously,  and  to  expedite 
computation,  the  following  notation  and  assumptions  are  introduced: 


1  - 
E 


a 


1 


l-o- 


1 

V 


a  =  1,  r  =  l  ,  b 


(92) 
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Equation  (92)  automatically  correctly  nondimensionalizes  all  the  stress 
expressions.  The  nine  equations  then  become 

2C2  +  C3  =  0 
2C5  +  6C7  +  4C8  =  0 


2C5  +  6C6  -  "6C7  -  2C8  =  ° 


6C_  4C 

2C  +  — -  +  — - 
5  4  2 

V  V 


rr 


2C0  +  -£•  =  r 
2  2  rr 

V  0 


2  6C7  2C8 

2C  +  6C  T]  ~  — —  ~  — =  r 
5  6  '  4  2  r0 
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TJ  0 
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-Trri  <!‘ai>  +  Tr01  <f  "  °1>J  ^ 


2C  2(1  -  <r  )C 
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„)  +  fr  (-  -  2a  )  -  t  (~  ”  2a  )|  ’ 
V  y  rr1  v  3  V  r^  v3  l'J  ^ 


(93) 
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Equation  (91)  is  solved  for  C.  . . .  C„,  t  ,  t  ,  and  r  ,  with 
'  2  8  rrQ  rr^  r 

r\  =  1.05,  1.1,  1.2  and  two  cases  of  material  properties  corresponding  to 

the  soft  and  stiff  cylindrical  lining  of  the  dynamic  analysis. 

Soft  Cylindrical  Lining 


v  =  2. 9,  cr  =  .25,  a  =  .20 

J.  Z 

Stiff  Cylindrical  Lining 


v  =  .  31,  cr  =  .25,  cr  =  .30 

X.  Li 


The  values  of  constants  C.  . . .  C.,  r  ,  t  ,  and  r  for  these 

2  8  rr^  rr,  r0, 

0  1  1 

two  cases  are  tabulated  in  Table  I. 


Table  I 


Case  I  v  =  2.9  cr  =  .25  cr  =  .20 

1  Li 

V 

C2 

°3 

C5 

C6 

C7 

°8 

T 

rro 

T 

rri 

T 

rei 

1.05 

-.112 

.224 

.053 

.002 

.057 

-.112 

-.021 

.019 

.041 

1.1 

-.116 

.232 

.053 

.003 

.060 

-.116 

-.  040 

.033 

.077 

1.2 

-.123 

.246 

.055 

.004 

.063 

-.122 

-.075 

.047 

.133 

Case  II  v  =  .  31  a  =  .25  a2  =  •  30 

V 

C2 

°3 

C5 

C6 

C7 

°8 

r 

rro 

r 

rri 

T 

r*l 

1.05 

-.1004 

2.008 

.596 

-.052 

.491 

-1.035 

-.187 

.138 

.298 

1.1 

-.903 

1.  806 

.595 

-.069 

.458 

-  .984 

-.313 

.187 

.442 

1.2 

-.775 

1.551 

.572 

-.071 

.429 

-  .930 

-.474 

.197 

.578 
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Using  the  values  of  the  constants  tabulated  in  Table  I,  it  is  possible  to 

evaluate  the  stresses  and  displacements  at  the  boundaries  of  the  stiff  and  soft 

cylinder.  For  purposes  of  this  report,  r*  ,  r*  ,  and  r*  are 

(2)  (1)  rr(l) 

r  -  a  r  =  b  r  =  b 

calculated  for  each  case.  The  (*)  notation  is  consistant  with  the  nondimensional 
notation  used  in  the  dynamic  analysis. 


T*ee  =  2C2  "  c3  +  (2C5  +  12C6  +  6C7>  cos  20 

(2) 

r  =  a 

r*  =  t  1  +  r  =  — - —  -  1  +  (2  -  4cr  )  cos  2  0  -  r 

00  00  00  1  -  a  v  1;  rr 

(1)  ,  (1)  (1)  1  L  J 

r  =  b  r  =  b  r  =  b 


T* 


rr 


(1) 


r  =  b 


T  + 

rro 


cos  2  0 


( T  -2  T  )  COS  2  0 
rr1  r0  7 


(94) 


These  values  for  0  =  7r,  tt/2  and  tj  =  1.05,  1.1,  1.2  are  plotted 
on  Figs.  3,  4,  5,  6,  8,  9,  10  and  11  at  a^a  =  0;  this  represents  the  infinite 
wave  length  or  static  solution. 
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